;  I  L  D 


AD-A230  748 


TATION  PAGE 


form  Approved 
0M8  No.  0704-0  r 83 


to  t  fvou#  pm  tmomym.  inctuOinQ  tf>«  tim«  for  inctrua*ort«.  M«fcn««g  axiviing  4«ta  tourcot. 

rwmg  mo  ooiloction  of  tnform«t*oo  S«Ad  commonto  rogordmg  mis  bordsn  s«i«nsts  or  sny  othor  ssooct  of  (f>« 
«sn.  to  Wssnngton  Hssoqwsnsrs  Sookss.  Otrsctorsto  for  Irrformstion  Oporsttons  ■f>0  Asoom.  1215  Jofforson 
•  of  Msns^stnsnt  sno  Sudgot.  Psoor-worli  Rsduaion  Profod  <0704-01 88l.  Wtshmgton.  OC  20903 


■  oianK/  2.  REPORT  DATE 

November  1990 


4  TITLE  AND  SUBTITLE 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final  from  10/1/89  Tn  )0/91/q 

I  5  FUNDING  NUMBERS 


Heat  rransfer  Predictions  of  Hypersonic  Transitional 
F  lows . 

6  AUTHOR(S) 

Joelle  M.  Champnev 


3x)- 

C--  ooo¥- 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Applied  d  Theoretical  Mechanics.  Inc.  (ATM) 
-501  Sequoyah  Road 
Oakland,  CA  94605 


9.  SPONSORING/MONITCRING  AGENCY  NAME(S)  AND  AD0RESS(ES) 

Air  Force  Office  of  Scientific  Research 
Bollint;  Air  Force  Base,  D.C.  20332 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


ATM-TR-90-025 


10.  SPONSORING/MOWTORING 
AGENCY  REPORT  NUMBER 

Xl>0  7//hl 


11  SUPPLEMENTARY  NOTES 


1 2a.  0ISTRIBUTI0N/AVAILA8IUTY  STATEMENT 


Di.stribuC  ion  unlimited 


JDi'Sm'TUTIOK  STfllTEMprr  A 
,'Af  fpr  p'-tbSe 


l^g^1 


JAN09  199I 


D 


13  abstract  (Maximum  200  'tvordsi 

.5  numerical  model  to  solve  transition  process  observed  in  hypersonic  flows 
over  cones  has  been  developed.  Low  Reynolds  number  two-equation  turbulence 
models  were  employed  with  a  production  term  modification  (PTM)  technique. 
The  approach  determined  the  extent  of  the  transition  zone.  The  onset  of 
transition  was  imposed  using  experimental  measurements  when  available.  l-Hien 
not  available,  the  onset  of  transition  was  determined  by  a  stability 
criterion  which  is  related  to  the  Bushnel 1-Reshotko  transition  criterion. 
The  PTM  technique  was  incorporated  into  a  NASA-Ames  implicit  Reynolds 
.iveraged  Navicr  Stokes  solver,  the  TURF  code,  and  tested  for  transitional 
livpersonii:  flows  over  flat  plates.  The  model  parameters  were  tuned  as  a 
] unction  of  free  stream  Mach  number.  The  PTM  technique  was  also  tested  for 
transitional  hypersonic  flows  over  sharp  cones  and  blunt  cones  for  a 
'.  irioLv  of  flow  conditions.  Comparisons  of  computed  heat  transfer  with 
exper  iment.i  1  measurements  are  shown  to  be  good.  ^  '' 


14  SUBJECT  TERMS 


i L  ion- hvpe r son ic -numeric a  1-  imp !  ic i t-two-equa t ion  model 


15  NUMBER  OF  PAGES 


16  PRICE  CODE 


17  SECURITY  CLASSIFICATION  18  SECURITY  CLASSIFICATION  19  SECURITY  CLASSIFICATION  20  LIMITATION  OF 

OF  report  of  this  page  of  ABSTRACT  ABSTRACT 


. :  0  '  1  .  i  s ,  ^  i  i  i  t.  n 


N5N  7540  01  280  5500 


nni  1  .i.T.s  i  i  1  (1 


nnc  1  .iss  i  t  ied 


i  !  1  .  '  ^ 


nn ] imiied 

Standard  Form  298  |R«w  2  09) 

f*r««criCMd  W  4NS<  Si*  2  39  1  0 
296  102 


HEAT  TRANSFER  PREDICTIONS  OF  HYPERSONIC  TRANSITIONAL  FLOWS. 


J.M.  Champney 

Applied  &  Theoretical  Mechanics,  Inc.  (ATM) 
A501  Sequoyah  Road 
Oakland,  CA  9A605 


Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
(AFSR) ,  under  Contract  F49620-90-C-0004 .  The  United  States 
Government  is  authorized  to  reproduce  and  distribute  reprints  for 
governmental  purposes  notwithstanding  any  copyright  notation 
hereon. 


The  views  and  conclusions  contained  in  this  document  are  those  of 
the  author  and  should  not  be  interpreted  as  necessarily 
representing  the  official  policies  or  endorsements,  either 
expressed  or  implied,  of  the  Air  Force  office  of  Scientific 
Research  or  the  U.S.  Government. 


TABLE  OF  CONTENTS 

SECTION  Page 

TABLE  OF  CONTENTS  2 

ILLUSTRATIONS  3 

FOREWORD  5 

ABSTRACT  6 

1 .  INTRODUCTION  7 

2.  TECHNICAL  APPROACH  8 

2.1  The  numerical  model.  8 

2.2  The  PTM  technique.  9 

2.3  Conversion  of  noise  level  into  free  stream  turbulence.  14 

2.4  The  stability  criterion.  15 

2.5  The  law  of  the  wall.  15 

3.  RESULTS.  16 

3.1  Flat  plates.  16 

3.1.1  Incompressible  regime.  16 

3.1.2  Richards  flow.  17 

3.1.3  Keener  flow.  18 

3.2  Sharp  cones.  19 

3.2.1  Martellucci  experiment,  sharp  cone.  19 

3.2.2  Muir  and  Trujillo  experiment,  sharp  cone.  19 

3.3  Blunt  cones.  20 

3.3.1  Martellucci  experiment,  blunt  cone.  21 

3.3.2  Muir  and  Trujillo  experiment,  blunt  cone.  21 

4.  CONCLUSIONS.  22 

5.  REFERENCES.  23 


2 


ILLUSTRATIONS 

FIGURE  NO.  Page 

1.  Flow  disturbances  in  Wind  Tunnels.  25 

2 .  Variation  of  "A"  with  free-stream  turbulence 

intensity  for  the  Jones-Launder  model.  26 

3.  Variation  of  "B”  with  free-stream  turbulence 

intensity  for  the  Jones-Launder  model.  26 

4.  Free-stream  rms  pressure  fluctuations  divided  by 

mean  static  pressure  versus  Mach  number.  27 

5.  Skin  friction  versus  Reynolds  number  based  on  streamwise 

distance,  computed  by  Schmidt  and  Patankar  with  the 
Jones-Launder  model  without  the  PTM  technique.  28 

6.  TURF  computed  skin  friction  versus  Reynolds  number 

based  on  streamwise  distance,  with  the  Jones-Launder 
model  and  without  the  PTM  technique.  28 

7 .  Skin  friction  results  of  Schmidt  and  Patankar 

obtained  with  the  Jones-Launder  model  and  the  PTM 
technique.  29 

8.  TURF  computed  skin  friction  with  the  Jones-Launder 

model  and  the  PTM  technique.  29 

9.  Profiles  of  turbulent  kinetic  energy  for  a  flat  plate 

flow  at  Mach  8.2  -  No  PTM  used.  30 

10.  Profiles  of  turbulent  kinetic  energy  for  a  flat  plate 

flow  at  Mach  8.2  -  PTM  used.  30 

11.  Velocity  profiles  along  flat  plate,  M=8.2  -  No  PTM.  31 

12.  Velocity  profiles  along  flat  plate,  M=8.2,  PTM  used.  31 

13.  Comparisons  between  computed  and  experimental  heat 

transfer.  Re/in  =  560,000.  32 

14 .  Comparisons  between  computed  and  experimental  heat 

transfer,  M=8 . 2 ,  showing  TUE  effect.  33 

15.  Comparisons  between  computed  and  experimental  heat 

transfer.  34 

16.  Comparisons  of  computed  and  measured  Stanton  number 

for  flow  over  a  sharp  cone  at  Mach  8.  35 

17.  Streamwise  variation  of  Stanton  number,  no  transition 

model  used.  36 


3 


37 


18.  Comparisons  between  computed  Stanton  number  and 

experimental  measurements.  The  PTM  technique  is  used. 


19.  Streamwise  variation  of  Stanton  number,  the 

intermittency  factor  approach  of  Adams  is  used.  38 

20.  Variation  of  Reynolds  number  based  on  momentum 

thickness  versus  Rex.  39 

21.  Blunt  cone  mesh  in  the  nose  region.  40 

22.  Density  contours  for  the  blunt  cone  at  Mach  8.  40 

23.  Entropy  contours,  blunt  cone  at  Mach  8.  41 

24.  Enlargement  of  entropy  contours,  blunt  cone  at 

Mach  8.  41 

25.  Calibration  of  A  and  B  in  the  PTM  model.  42 

26.  Variation  of  Stanton  Number  with  free-stream 
Reynolds  number  for  a  blunt  cone,  RN=0.1  in.  - 

no  transition  model  used.  43 

27.  Variation  of  Stanton  Number  with  free-stream 
Reynolds  number  for  a  blunt  cone,  RN=0.1  in.  - 

PTM  model  used.  44 


4 


FOREWORD 


The  author  wishes  to  thank  Dr.  T.J.  Coakley  for  many  useful 
discussions  and  technical  advice.  The  author  is  grateful  to  the 
Applied  Computational  Fluid  Branch  at  NASA-Ames  Research  Center 
for  providing  the  computer  time. 


5 


ABSTRACT 


A  numerical  model  to  solve  transition  process  observed  in 
hypersonic  flows  over  cones  has  been  developed.  Low  Reynolds 
number  two-eqfuation  turbulence  models  were  employed  with  a 
production  term  modification  (PTM)  technique.  The  approach 
determined  the  extent  of  the  transition  zone.  The  onset  of 
transition  was  imposed  using  experimental  measurements  when 
available.  When  not  available,  the  onset  of  transition  was 
determined  by  a  stability  criterion  which  is  related  to  the 
Bushnell-Reshotko  transition  criterion.  The  PTM  technique  was 
incorporated  into  a  NASA-Ames  implicit  Reynolds  averaged  Navier 
Stokes  solver,  the  TURF  code  and  tested  for  transitional 
hypersonic  flows  over  flat  plates.  The  model  parameters  were 
tuned  as  a  function  of  free  stream  Mach  number.  The  PTM  technique 
was  also  tested  for  transitional  hypersonic  flows  over  sharp 
cones  and  blunt  cones  for  a  variety  of  flow  conditions. 
Comparisons  of  computed  heat  transfer  with  experimental 
measurements  are  shown  to  be  good. 
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1 .  INTRODUCTION 


This  research  effort  is  to  develop  a  numerical  approach 
to  predict  heat  transfer  variation  in  the  transition  zone  as 
observed  on  hypersonic  aircraft.  The  accurate  prediction  of  heat 
transfer  on  the  vehicle  airframe  has  been  recognised  to  be 
critical  in  the  development  and  design  processes.  Peak  heat 
transfer  rates  and  duration  of  heating  determine  the  mass  of 
thermal  protection  systems  on  spacecraft.  Heating  rates  are 
strongly  influenced  by  the  boundary  layer  transition  to 
turbulence,  causing  extremes  in  heat  transfer  and  skin  friction. 
Transition  is  the  process  by  which  a  laminar  boundary  layer 
becomes  turbulent.  This  process  is  complex  because  it  is 
influenced  by  a  variety  of  factors.  Reference  1  gives  a  review  of 
transition  theory  and  describes  various  techniques  employed  to 
analyse  transition. 

The  first  theoretical  idea  for  the  analysis  of  transition  is 
due  to  O.  Reynolds  and  Lord  Rayleigh  in  the  1880.  Fifty  years 
later,  theoretical  investigations  lead  to  the  stability  theory 
formulated  by  Prandtl's  school  in  about  1930.  The  next  forty 
years  saw  theoretical  and  experimental  investigations  to  examine 
transition  processes.  Over  this  period  and  up  to  the  present, 
summaries  and  review  papers  on  transition  appeared  but  hypersonic 
transition  received  little  theoretical  investigations,  as  a 
result  of  the  reduced  interest  for  hypersonic  flow  regimes.  Even 
subsonic/transonic  transition  received  reduced  attention  over  the 
last  fifteen  years  (ref.  2) .  Flight  test  data  and  results  of 
current  wind  tunnel  experiments  indicate  that  boundary  layer 
transition  has  a  major  influence  on  the  aerodynamic  behaviour  of 
hypersonic  vehicles  and  specifically,  has  a  major  effect  on; 

skin  friction  drag, 
surface  heat-transfer  rates, 

flow  separation  and  surface  pressure  distributions, 

extent  of  shock  boundary-layer  interaction, 

flow  separation  and  control  effectiveness, 

unsteady  flow  phenomena, 

vehicle  stability  and  control, 

surface  pressure  fluctuations, 

structural  fatigue, 

acoustic  noise, 

and  others  (see  Pate,  ref. 3). 

Since  the  work  of  Kovasnay  in  1953  (ref.  4) ,  the  response  of 
the  boundary  layer  to  environmental  disturbances  is  now  believed 
to  be  the  dominant  criterion  of  the  transition  process  (see  work 
by  Stetson,  ref.  5) .  In  flight,  disturbances  have  much  lower 
freestream  disturbance  levels  than  conventional  supersonic- 
hypersonic  wind  tunnels.  Experimentalists  in  hypersonic 
transition  have  devoted  much  effort  to  reduce  facility  free¬ 
stream  disturbances  to  low  levels,  approaching  those  in  flight 
(ref.  6).  The  paper  by  Dougherty  and  Fisher  (ref.  7)  reviews  wind 
tunnel  flight  data  correlations  and  is  a  useful  source  of  data. 
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Free-stream  disturbance  effects  need  to  be  incorporated  into 
a  transition  model.  The  present  work  is  an  attempt  to  model  the 
interaction  of  the  free-stream  disturbances  and  turbulence  with 
boundary  layers  on  flat  plates  and  cones,  making  use  of  experimental 
data  and  recent  advances  in  Computational  Fluid  Dynamics. 

Kovasnay  identifies  three  primary  disturbance  sources  in 
wind  tunnels:  vorticity  fluctuations  (free-stream  turbulence) , 
entropy  fluctuations  and  sound  waves.  The  paper  by  Pate  (ref.  3) 
presents  a  review  of  these  effects.  Figure  1  extracted  from  Pate 
paper,  summarizes  the  dominant  effects  of  flow  disturbances  as  a 
function  of  Mach  numbers  and  as  observed  in  wind  tunnels.  It  can 
be  seen  in  Figure  1  that  the  radiated  noise  is  the  dominant  free- 
stream  effect  of  the  boundary  layer  transition  processes  for  the 
Mach  numbers  of  interest  in  this  effort,  i.e.  greater  than  3.  The 
importance  of  temperature  fluctuations  is  not  well  understood  at 
the  present  time. 

In  this  research  effort,  the  transition  model  developed  by 
Schmidt  and  Patankar  (ref. 8)  for  incompressible  flows,  was  tested 
and  modified  for  transitional  hypersonic  flows  using  experimental 
data  with  the  most  influencial  factor  in  the  model  taken  to  be 
the  level  of  noise  converted  into  free-stream  turbulence  intensity. 
The  approach  employs  a  production  term  modification  (PTM)  which 
consists  of  a  modification  to  low  Reynolds  number  (LRN)  two- 
equation  models.  A  modification  to  the  production  term  in  the 
modeled  turbulence  field  equations  is  correlated  to  the  free- 
stream  turbulence  level  using  experimental  data.  The  modification 
does  not  affect  the  fully  turbulent  calculations  and  greatly 
improves  the  transition  predictions.  The  approach  is  discussed  in 
the  following  section. 


2.  TECHNICAL  APPROACH 
2.1  The  numerical  model. 

The  computer  code  selected  for  this  effort  is  the  NASA-Ames 
TURF  code  developed  by  T.J.  Coakley  (refs.  9  and  10) .  The  TURF 
code  is  an  implicit  finite  difference  code  for  solving  the 
Reynolds  averaged  compressible  Navier-Stokes  equations  and 
incorporates  several  zero-  and  two-equation  turbulence  models. 
The  code  solves  the  two-dimensional  (2D)  and  axisymmetric 
equations  using  a  finite  volume  approach  and  an  approximately 
factored  alternating  direction  implicit  (ADI)  algorithm.  Second- 
order  upwind  differencing  is  used  for  inviscid  flux  differencing 
(closely  related  to  the  flux  difference  splitting  approach  of 
Roe)  and  central  differencing  is  used  for  the  viscous  terms.  The 
implicit  operator  utilises  a  non-conservative  diagonal  form, 
first  order  upwind  differencing  of  inviscid  terms  and  an 
approximate  (diagonalised)  second-order  differencing  of  viscous 
terms . 
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ATM  examines  the  possibility  of  replacing  the  left  hand  side 
(LHS)  operator  of  the  TURF  code  by  a  lower-upper  symmetric 
successive  overrelaxation  (LU-SSOR)  technique,  as  developed  by 
Yoon  and  Jameson  (ref.  11) .  In  this  technique,  the  Navier-Stokes 
equations  are  differenced  and  not  factored  as  in  the  ADI  methods. 
The  technique  was,  at  first,  attractive  because  it  is  a  two 
factor  (or  sweep)  scheme  even  in  three  dimensions  (3D) ,  and  it  is 
theoretically  unconditionally  stable  in  3D.  This  was  found  not  to 
be  verified  for  simple  two  dimensional  flows,  it  was  not  possible 
to  use  a  large  time  step.  Various  forms  of  the  implicit  operator 
were  tested  using  at  first  u,  u+c,  u-c,  v,  v+c,  v-c,  for  the 
inviscid  eigenvalues  of  the  matrices  that  appear  in  the 
differencing  of  the  LHS  operator,  u  is  the  streamwise  (or  x) 
velocity  component,  v  is  the  normal  (or  y)  velocity  component,  c 
is  the  sound  speed.  The  computation  was  unstable  and  could  use 
only  a  small  time  step  corresponding  to  a  Courant-Friedrichs-Levy 
number  of  1.  If  each  eigenvalue  is  replaced  by  its  largest  value, 
then  the  time  step  can  become  infinite,  however,  the  rate  of 
convergence  to  a  steady  state  solution  is  slow,  slower  than  the 
original  TURF  upwind  ADI  algorithm.  T.J.  Coakley  and  George  Hwang 
at  NASA-Ames  Research  Center  have  pursued  extensive 
investigations  of  the  LU-SSOR  algorithm  that  lead  to  the 
implementation  in  the  TURF  code  of  a  Gauss  Seidel  line  relaxation 
technique  for  the  LHS  operator.  In  their  applications,  the  new 
algorithm  allows  an  infinite  time  step,  but  even  with  the 
infinite  time-step,  the  code  still  converges  in  roughly  the  same 
number  of  time  steps  than  the  original  TURF  code.  In  addition, 
the  Gauss  Seidel  technique  is  not  vectorizable  in  2D  on  the  NAS 
CRAY,  and  uses  more  computer  time  than  the  TURF  algorithm.  For 
theses  reasons  and  following  T.J.  Coakley  recommendations,  it  was 
decided  to  keep  in  the  TURF  code,  the  original  LHS  implicit  ADI 
operator. 


2.2  The  PTM  technique. 

The  purpose  of  this  section  is  to  describe  the  mathematical 
implementation  of  the  PTM  technique  in  the  low  Reynolds  number 
Jones-Launder  k-e  turbulence  model,  k  is  the  turbulent  kinetic 
energy  and  e  is  the  dissipation  rate  of  k. 

The  basic  relation  defining  the  turbulent  viscosity  is 


^T 


P  C 
^  P 
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p  is  the  density,  C^  is  a  constant  equal  to  0.09,  f^  is  a  low 
Reynolds  number  function  or  damping  function. 
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exp 


-3.4 


(1+0. 2R^) 


2 

is  the  turbulent  Reynolds  number  and  is  equal  to  k  /ve, 
y^is  the  kinematic  viscosity  equal  to  fi/p,  p  is  the  molecular 
viscosity  given  by  Sutherland  law. 


The  unsteady  transport  equations  for  k  and  e  are,  in  2D: 


|£k  ,  |_  {p.v-  .  P,- 


iei  +  L. 

at 


Pj^  is  the  production  of  k  and  is  equal  to:  -  (2/3)  pkD. 


au  ^  av 
^  ^k,k  ax  ay 


2  2  2 

S  =  (U.  .  +  U.  .)U.  .  -  4 

1,3  1,3^  1,3  3 

2  V  /VTc\  ^ 


f^  =  1.  ,  f^  =  l.-0.3exp(-l4)  ,  D  =  (g)  ,  E  =  2p^v  (|H) 


=  1.44,  =  1.92,  =  1.0,  =  1.3. 

t  is  the  time,  the  indices  i,  j,  k  indicate  the  i,  j,  k 
directions  and  can  be  equal  to  1  (i.e.,  x)  or  2  (i.e.,  y) . 

ATM  experimented  with  the  Chien,  Lam-Bremhorst  and  Jones- 
Launder  models  (refs.  12-14)  for  flat  plate  flows  and  did  not 
find  any  difference  between  the  models  that  could  allow  to  define 
one  of  the  models  as  the  "best”  model  for  this  study.  Consequently, 
in  the  applications  reported  here,  the  Jones-Launder  model  was 
used  because  it  is  the  oldest  two-equation  model  and  has  been 
widely  used. 

A  basic  difference  between  the  three  models  is  that  they 
each  use  different  LRN  functions  for  f  ,  f-,  f ^ ,  E,  D.  Another 
difference  between  the  three  models  is^that  in'^the  Chien  and 
Jones-Launder  models,  the  variable  e  is  a  modified  dissipation 
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function,  so  that  its  value  at  the  wall  can  be  set  to  zero.  In 
the  Lam-Bremhorst  model,  e  retains  its  original  definition  which 
leads  to  zero-gradient  boundary  condition  on  e  at  the  wall. 

Initial  and  far  field  conditions.  The  stability  of  a  two-equation 
model  computations  is  usually  related  to  the  free  stream  values  of 
the  turbulent  variables:  k^,  and  ,  »  denotes  free  stream 
values.  The  following  values  are  used'^For  initial  conditions: 


k  -  i  (TUE*U  )^  ,  =  0.01  fi  ,  e  =  RSR*Mt,/(C  pk^) 

TUE  is  the  free  stream  turbulence  level  that  can  vary  from  0.1% 
to  6%.  The  stability  of  the  computation  is  usually  related  to 
the  value  of  e  .  With  the  Jones-Launder  model,  if  e  is  too 
large,  turbulence  is  destroyed  and  the  flow  remains“laminar ,  if 
e  is  too  small,  the  computation  becomes  unstable.  There  is 
therefore,  an  optimum  value  of  e  .  In  our  computations,  it  was 
found  that  an  appropriate  value  of  e  corresponds  to  the  constant 
RSR  equal  to  0.001  or  0.01.  The  k-e  solution  does  not  depend  upon 
these  values.  In  the  far  field,  k  and  e  are  set  to  their  initial 
values.  In  the  flow  field,  they  are  prevented  from  becoming 
smaller  than  k  and  e  . 

CO  CO 

Generally,  two-equation  models,  when  applied  to  transitional 
flows  over  flat  plates,  predict  laminar  flow  followed  by  fully 
turbulent  flow.  In  the  present  applications,  ATM  found  that  the 
two-equation  models  predict  fully  turbulent  flow  starting  at  the 
leading  edge  of  flat  plate  flows,  when  crude  initial  conditions 
are  used,  i.e.  the  flow  variables  are  set  to  constant  values.  If 
initial  conditions  are  refined  using  known  analytical  solutions 
for  the  flow  variables,  the  two-equation  models  will  predict 
laminar  flow  followed  by  fully  turbulent  flow.  In  all  cases,  the 
transition  zone  is  very  narrow.  This  is  translated  by  a  sharp 
increase  of  skin  friction  and  heat  transfer.  The  profiles  of 
turbulent  kinetic  energy,  k,  along  flat  plates,  in  the  laminar 
region  show  a  very  slow  increase  from  zero  at  the  wall  up  to  the 
boundary  layer  edge  value.  Downstream,  the  k-profiles  undergo  the 
rapid  development  of  a  peak  inside  the  boundary  layer:  turbulent 
kinetic  energy  from  the  free  stream  is  convected  and  diffused 
into  the  boundary  layer.  In  this  process  the  production  of  k  in 
the  two-equation  model  becomes  significant  and  larger  than  the 
dissipation  of  k.  The  process  is  self  generating  causing  a  rapid 
increase  of  k.  It  is  necessary  therefore  to  modify  the  production 
term  of  k  in  the  k  and  e  equations  to  improve  predictions  of 
transition  processes. 

The  PTM  technique  is  employed  in  conjunction  with  a  two- 
equation  model.  The  PTM  technique  limits  the  growth  rate  of  P^^, 
using  a  linear  relationship  depending  on  Pj^: 
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A  +  B 
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A  and  B  have  the  dimension  of  the  inverse  of  time. 


Another  way  to  compute  the  growth  rate  of  Pj^  is  to  use  the 
following  equation: 
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where  At  is  the  computational  time  step  (=dx/u  where  dx  is  the 
stepsize  in  the  streamwise  direction) . 

P?^  is  the  production  of  k  at  time  n  resulting  from  the  use 
of  the^PTM  technique. 

pj^^^  is  obtained  from  the  product  of  the  Reynolds  stresses 
with  the  gradients  of  the  mean  velocity. 

The  PTM  technique  determines  the  value  of  P.  at  x+dx  to  be 
used  over  the  next  step  in  the  solution,  by  adding  to  the  value 
of  P.  at  X,  the  smallest  increment  (AP.  )  for  the  growth  rate  of 
Pj^,  obtained  from  equation  1  or  2. 

Equation  1  is  employed  in  the  transition  zone,  characterised 
as  the  zone  where  the  calculated  Reynolds  number  based  on 
momentum  thickness,  R^a,  is  greater  than  a  critical  Reynolds 
number,  R^^  ^  (equal  to  162  for  incompressible  flows) . 

If  R  is  less  than  R  „  ,  the  growth  rate  of  P  is  set  to 

zero  and  ?ne  two-erruation  model  will  predict  laminar^flow.  When 
the  growth  rate  of  P.  calculated  by  (1)  becomes  larger  than  the 
growth  rate  of  Pj^  given  by  Equation  2,  Equation  1  is  not  used 
anymore,  leaving^the  two-equation  model  unmodified  by  the 
transition  model,  in  fully  turbulent  flow. 

In  Equation  1,  the  two  ijarameters,  A  and  B,  are  made 
dependent  upon  the  free-stream  turbulence  intensity  and  were 
determined  using  the  correlation  of  Abu-Ghannam  and  Shaw  (ref. 

15) ,  which  is  a  criterion  for  R  'o  determine  the  start  and  the 
end  of  transition.  ® 

The  A  and  B  parameters  were  derived  in  reference  8  and  are 
given  below. 
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Leip=100*TuE,  B=BMO^^  and^  =  A»10^ 


12.7536*p^ 
O.lZlS^p'^ 

6.0  <p  loge(B)  =  1.950  +  0.1573*(j) 


0.0  <  9  <  2.0 


logeiB)  =-5.8084  +  2.995*p 


2.0  <  p  <  6.0 


B  =  18.738  -  26.8085*(})  + 
-  2.1I52*(!)^  + 


0.0  <  (!)  <  6.0 


A  =  12.2266  -  1.7904  *())  -  2.4229*p“ 
+  0.57595*(!)^  -  0.0365*p'^ 

A  =-7.5-0.19*(!) 


A  and  B  are  the  non  dimensional ized  A  and  B  parameters  with 
respect  to  local  free  stream  conditions  involving  the  density, 
the  velocity  and  the  molecular  viscosity: 


Am, 


A  = 


P 

'^e  e 


B  M 


B  = 


3  tt6 


Figures  2  and  3  give  the  variation  of  A  and  B  as  a  function  of 
TUE  for  the  Jones-Launder  model. 

From  Equations  1  and  2,  it  can  be  noted  that  the  time 
dependency  is  built  in  the  model  to  slow  down  the  transition 
process.  The  time  scale  is  related  to  the  local  velocity.  The  PTM 
technique  verify  some  experimental  evidences  (ref.  8): 

-  Free  stream  turbulence  influences  the  transitional  flow. 

-  When  the  free-stream  turbulence  intensity  increases,  the 
transition  region  moves  upstream,  skin  friction  and  Stanton 
number  increase  too. 

-  In  accelerating  free  stream  flows  (favorable  pressure  gradient) 
the  transition  moves  downstream,  the  transition  zone  is  longer 
because  of  the  stabilizing  effects  of  accelerated  freestream  flows. 

-  The  process  of  transition,  once  started,  evolves  at  a  finite 
rate,  the  rate  at  which  Pj^  changes  is  limited. 
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Schmidt  and  Patankar  compared  their  incompressible  model 
against  experimental  measurements  for  34  low  speed  flow 
experiments.  In  most  cases,  the  agreement  of  skin  friction  and 
heat  transfer  is  good,  if  not  excellent.  In  the  present  approach, 
compressibility  effects  are  treated  in  terms  of  Favre  (mass- 
averaged)  variables  and  the  gas  is  assumed  perfect  with  a 
constant  "gamma",  therefore,  Favre  averaging  is  appropriate.  The 
use  of  Favre  averaging  is  a  controversial  issue  in  the  turbulence 
modeling  of  hypersonic  flows  because  hypersonic  flows  involve 
multiple  reacting  species  interacting  with  turbulence.  Favre 
averaging  tends  to  oversimplify  the  partial  differential 
equations,  important  terms  might  disappear.  However,  this  very 
complex  issue  is  not  addressed  in  this  work. 


2.3  Conversion  of  noise  level  into  free  stream  turbulence. 

To  apply  two-equation  turbulence  models,  a  value  of  the 
free-stream  turbulence  level,  is  needed  in  the  far  field  to  start 
the  computation  of  the  turbulent  kinetic  energy  k  and  to  set  the 
far  field  boundary  conditions.  In  addition,  the  PTM  model 
includes  two  constants  that  depend  upon  the  free  stream 
turbulence  level.  In  hypersonic  flow,  the  noise  level  is  assumed 
to  be  the  dominant  influencial  factor  and  is  related  to  the  free 
stream  fluctuating  pressure.  The  free  stream  turbulence  is 
related  to  the  free  stream  velocity  fluctuation.  The  noise  level 
is  now  measured  by  experimentalists  in  wind  tunnels  and  it  is 
therefore  desirable  to  relate  the  noise  level  to  a  free  stream 
turbulence  level  to  perform  numerical  simulations  with  two 
equation  models.  Laufer  (ref.  16)  derives  a  relation  assuming 
that  the  free  stream  disturbance  is  a  plane  wave  propagating  at  a 
Mach  number  relative  to  the  free  stream.  Using  the  unsteady 
Bernoulli  equation,  an  equation  between  fluctuating  velocity  and 
fluctuating  pressure  is  obtained.  This  relation,  used  by  Shamroth 
and  MacDonald  (ref.  17) ,  involves  a  factor  treated  as  an  integral 
of  a  space-time  correlation  function  which  is  difficult  to 
evaluate  and  crude  assumptions  have  to  be  made.  Therefore,  a 
simple  analysis  to  relate  the  fluctuating  pressure  and  TUE  was 
performed  using  a  one  dimensional  perturbation  technique  of  the 
steady  Bernoulli  equation.  A  perfect  gas  is  assumed,  temperature 
fluctuations  are  neglected.  The  steady  inviscid  Bernoulli 
equation  is  assumed  verified  in  wind  tunnels  in  the  free  stream 
and  is  written  as: 


E 

P 


+ 


1 

2 


constant 


The  variables  pressure,  p,  density,  p,  and  magnitude  of  the 
velocity,  q,  are  perturbed  from  their  free  stream  values: 

P  =  P  +  P'  /  P  =  +  P'  ,  =  (U  +u')^  +  v'^  +  w'^ 

OO  00  00 

with  u'  >>  V'  and  u'  >>  w' .  The  perturbed  Bernoulli  equation 
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becomes : 
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With  this  crude  analysis,  the  experimental  data  of  rms 
pressure  fluctuations  versus  Mach  number  as  given  in  Figure  4  can 
be  used  to  obtain  the  fluctuating  pressure  for  a  given  Mach 
number  and  subsequently  the  value  TUE  needed  in  the  numerical 
simulations . 


2.4  The  stability  criterion 

The  PTM  technique  employs  a  stability  criterion  stating  that 
for  a  Reynolds  number  based  on  momentum  thickness,  R  ,  less  than 
a  critical  value,  R^g  the  production  term  of  k  is^set  to  zero. 

For  incompressible  flows,  R  q  ^  was  selected  to  be  125  (ref. 
8) .  For  hypersonic  flows,  the  sttbllity  criterion  is  not  known, 
at  the  present  time.  Reshotko  et  al.  (ref.  2)  recommend  for  2D 
and  axisymmetric  flows:  R  ^  =  K  a/M  =  150,  where  M  is  the 

local  edge  Mach  number.  ' 

This  criterion  is  used  to  determine  the  beginning  of 
transition,  i.e.,  the  location  of  minimum  heat  transfer  level.  In 
the  PTM  model,  the  growth  of  P.  can  lag  after  being  triggered  by 
the  stability  criterion,  therefore,  R  _  may  need  to  be  smaller 
than  150.  In  addition,  it  would  be  desifSble  to  include  in  the 
PTM  stability  criterion  the  variation  of  TUE.  This  variation  will 
be  accurately  included  when  more  hypersonic  transitional 
experimental  data  with  measured  noise  level  (to  be  converted  into 
TUE  level)  will  become  available.  For  the  present  time,  it  is 
recommended  to  use  R^  equal  to  125  in  the  PTM  technique.  It 
should  also  be  noted^that  the  quantity  R  _  is  not  easy  to  compute 
accurately  in  particular  for  hypersonic  flows  over  cones,  because 
the  edge  of  the  boundary  layer  and  therefore  the  momentum 
thickness,  8,  is  difficult  to  locate  for  these  flows. 


2.5  The  law  of  the  wall 

The  TURF  code  incorporates  the  capability  of  using  the  Law 
of  the  Wall  boundary  condition  (LOW  BC)  with  the  high  Reynolds 
number  form  of  two-equation  models.  This  type  of  boundary 
conditions  is  described  in  references  18  and  19.  Its  previous 
applications  for  fully  turbulent  flows  with  large  separation  (in 
particular,  step  flow)  were  quite  successful  regarding  numerical 
efficiency.  For  example,  without  LOW  BC,  the  step  flow  of  ref.  19 
could  not  be  computed  because  for  this  flow  the  use  of 
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integration  to  the  wall  boundary  condition  lead  to  prohibitively 
expensive  computer  runs.  In  the  LOW  BC,  the  wall  shear  stress  is 
determined  implicitly  assuming  that  the  computed  velocity  at  the 
first  cell  center  verifies  a  logarithmic  profile.  This  assumption 
allows  the  use  of  a  larger  spacing  in  the  direction  normal  to  the 
wall  corresponding  to  a  dimensionless  distance  y  at  least 
greater  than  10.0,  while  with  integration  to  the  wall  BC,  y  must 
remain  approximately  less  or  equal  to  1.0. 

Numerical  experimentations  were  made  to  investigate  the 
feasibility  of  using  LOW  BC  for  transitional  hypersonic  flows 
over  flat  plates.  The  results  of  heat  transfer  were  found 
sensitive  t^  the  wall  spacing.  For  a  dimensionless  distance  to 
the  wall,  y^,  of  about  10,  the  peak  of  heat  transfer  was  too 
high.  For  y^  of  about  100,  it  was  too  low.  Therefore,  by 
adjusting  y  ,  agreement  with  experiment  could  be  reached 1  The 
inconveniency  of  the  LOW  BC  is  to  use  a  velocity  profile  which  is 
logarithmic,  this  is  not  true  in  the  transition  region.  The 
behaviour  of  the  velocity  profile  in  the  transition  region  of 
high  speed  flows  is  not  well  known  at  the  present  time.  Much  more 
theoretical  work  is  needed  to  make  the  LOW  BC  applicable  to 
transitional  hypersonic  flows.  The  use  of  integration  to  the 
wall  with  LRN  two-equation  models  is  recommended  and  was  employed 
in  this  study. 

3.  RESULTS 

Numerical  simulations  of  hypersonic  transitional  flows  over 
flat  plates,  sharp  cones  and  blunt  cones  were  performed  with  the 
TURF  code  and  are  reported  in  this  section. 

3.1  Flat  plates. 

3.1.1  Incompressible  regime. 

The  implementation  of  the  PTM  technique  in  the  TURF  code  was 
tested  for  flows  over  flat  plates  in  the  incompressible  regime. 
Subsonic  boundary  conditions  are  used.  The  mesh  size  consist  of 
100  uniformly  spaced  points  in  the  streamwise  direction  and  40 
exponentially  stretched  points  in  the  normal  direction.  Figures  5 
and  6  compare  skin  friction  results  of  Schmidt-Patankar  with  the 
TURF  results,  using  the  Jones-Launder  model,  without  the  PTM 
technique  for  free  stream  turbulence  levels  ranging  from  1%  to  6%. 
Figures  7  and  8  compare  skin  friction  results  of  Schmidt- 
Patankar  with  the  TURF  results,  using  the  Jones-Launder  model  and 
the  PTM  technique.  Figure  6  shows  that  the  TURF  code  with  the 
Jones-Launder  model  without  the  PTM  technique  overpredicts  the 
level  of  skin  friction  in  the  turbulence  region  for  the  highest 
levels  of  free-stream  turbulence  (TUE  >  3%).  In  Figure  8,  the 
peak  of  skin  friction  is  smoothed  by  the  PTM  technique  resulting 
in  a  spreading  of  the  rise  in  skin  friction  and  therefore  of  the 
transition  zone.  TURF  results  of  Figures  6  and  8  agree  with  the 
Schmidt-Patankar  results  only  satisfactorily.  They  could  be 
improved  with  a  finer  mesh.  Such  computations  require  300  CPU 
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seconds  on  the  NAS  CRAY,  with  a  finer  mesh  computer  time  could  be 
multiplied  by  10  per  case  (i.o.,  per  TUE  value),  this  was  outside 
the  computer  ressources  of  this  project. 

3.1.2  Richards  flow. 

In  1960,  B.E.  Richards  (ref.  20)  conducted  an  experimental 
investigation  of  natural  transition  and  turbulent  boundary  layers 
on  a  cold  flat  plate  in  hypersonic  flow  in  the  Imperial  College 
Gun  Tunnel.  The  heat  transfer  equipment  as  described  by  Holden 
(ref.  21)  was  employed.  Numerical  simulations  of  this  experiment 
were  performed  with  the  following  test  conditions: 

angle  of  incidence=  0®, 

free  stream  Mach  number,  II  =  8.2, 

stagnation  temperature,  T^  =  775®K, 

free  stream  unit  Reynolds^number ,  u  /v  per  in.  =  5.6x10^, 
model  :  a  1  foot  long  flat  plate.  “  " 

ratio  of  wall  temperature  to  recovery  temperature,  T^/T^  =  0.44. 

The  mesh  consists  of  100  cells  uniformly  spaced  in  the  streamwise 
direction  and  40  cells  exponentially  stretched  in  the  normal 
direction.  The  first  spacing  in  the  normal  direction  is  a 
function  of  the  free  stream  Reynolds  number. 

The  boundary  conditions  are: 

-  no  slip  on  the  flat  plate  surface,  integration  to  the  wall  is 
applied.  The  wall  temperature  is  imposed  to  0.44*T  . 

-  inflow:  all  flow  variables  are  fixed.  ^ 

-  outflow:  all  flow  variables  are  extrapolated. 

Implicit  boundary  conditions  are  used  for  the  LHS  operator. 

Such  mesh  and  boundary  conditions  are  used  for  all  the  flow 
simulations  reported  here. 

The  beginning  of  transition  is  imposed  to  be  0.58  inches. 

Figures  9  and  10  compare  the  results  of  turbulent  kinetic 
energy,  TKE,  profiles  along  the  flat  plate  without  and  with  the 
PTM  technique.  In  Figure  9,  a  sudden  change  of  TKE  is  observed  at 
transition  onset,  then  the  TKE  profiles  grow  and  flatten  towards 
a  fully  turbulent  (narrower)  profiles.  In  Figure  10,  the  TKE 
profiles  grow  monotonically  in  the  the  transition  region  and 
converge  toward  a  fully  turbulent  profile.  The  same  qualitative 
behaviour  was  observed  by  Schmidt  and  Patankar  for  incompressible 
flows.  Figures  11  and  12  show  the  comparisons  of  velocity 
profiles  along  the  flat  plate  without  and  with  the  PTM  technique. 
The  profiles  behave  in  the  same  qualitative  way  than  the  TKE 
profiles . 

Figure  13  compare  the  results  of  TURF  heat  transfer  with 
Richards  experimental  measurements,  for  a  free  stream  turbulence 
level  of  1%.  In  this  Figure,  we  observe  that  when  the  transition 
model  is  not  used,  the  Jones-Launder  (J-L)  model  predicts 
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reasonably  well  the  level  of  heat  transfer  peak  but  the 
transition  zone  is  too  narrow.  When  the  J-L  model  with  the  PTM 
technique  is  used  with  A  and  B  values  as  derived  for  incompressible 
flows,  the  predicted  heat  transfer  is  id^^ntical  to  the  one 
predicted  by  the  J-L  model  without  the  PTM  technique.  Hovewer, 
the  PTM  technique  contributed  to  improve  the  rate  of  convergence 
of  the  computation.  Such  behaviour  of  the  PTM  solution  is  caused 
by  a  too  large  value  of  the  PTM  constant  B.  In  Figure  13,  results 
are  shown  for  A  and  B  divided  by  M.  In  this  case,  the  peak  of 
h^at  transfer  is  too  low.  When  A  is  not  modified  and  B  divided  by 
,  we  observed  a  delay  in  the  beginning  of  transition  but  the 
peak  of  heat  transfer  is  still  too  low.  In  Figure  13,  in  the 
laminar  regime,  the  computed  heat  transfer  is  too  low.  Various 
mesh  spacings  in  the  normal  and  streamwise  direction  were  tested 
without  any  effect  on  the  results.  This  too  low  level  of  heat 
transfer  needs  to  be  investigated,  it  could  be  caused  by 
experimental  or  numerical  errors,  or  to  a  physical  aspect  of  the 
flow  not  modeled  by  the  numerical  simulation. 

In  Figure  14,  the  effect  of  varying  TUE  (TUE=2%)  was 
investigated  for  the  cases  where  the  PTM  constants  A  and  B  are 
both  divided  by  M,  and  for  the  case  where  B  alone  is  divided  by 
.  In  the  first  case,  the  heat  transfer  has  increased  but  the 
transition  zone  is  narrower  than  for  the  case  with  TUE=1%,  the 
second  case  shows  some  improvement.  In  all  cases,  the  PTM 
technique  does  not  predict  the  experimentally  observed  slow  rise 
of  heat  transfer  starting  upstream  in  the  laminar  regime. 

3.1.3  Keener  flow. 

Starting  in  the  early  sixties,  Earl  Keener  and  Thomas  Polek 
(ref.  22)  performed  a  series  of  wind  tunnel  experiments  at  Ames 
Research  Center  to  measure  heat  transfer  on  flat  plates  for 
hypersonic  flows.  The  results  although  not  published  were  made 
available  for  this  study. 

In  this  experiment,  air  is  heated  to  temperatures  ranging 
about  670  to  1170  degrees  Kelvins.  This  is  necessary  to  prevent 
air  condensation  that  will  occur  otherwise,  causing  a  transitional 
two-phase  flow.  The  test  section  pressure  is  about  13  pounds  per 
square  feet.  The  gas  is  assumed  perfect.  The  molecular  viscosity 
is  determined  as  a  function  of  temperature  by  the  Sutherland  law. 

The  flow  conditions  are: 

free  stream  Mach  number,  M  =  7.4, 

free  stream  temperature  =  17 3 «  Rankines, 

wall  temperature  =  547®  Rankines,  corresponding  to  a  ratio  of  wall 
to  adiabatic  wall  temperature  of  0.3. 

The  flat  plate  length  is  1  meter.  TUE  was  set  to  1%. 

Figure  15  shows  results  of  heat  transfer  variation  along  the 
flat  plate.  The  transition  onset  is  imposed  and  equal  to  50  cm. 

In  the  fully  laminar  case  for  a  free  stream  Reynolds  number.  Re, 
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per  meter  of  1,890,000,  agreement  between  experimental  and 
computation  is  excellent,  indicating  that  the  computation  is 
numerically  accurate  regarding  the  mesh  size  and  assumptions 
(Sutherland  law,  perfect  gas  ...)•  For  higher  Reynolds  numbers,  in 
the  laminar  regime,  the  computed  heat  transfer  is  too  low.  The  J- 
L  model  without  and  with  the  PTM  technique  (A  and  B  equal  to 
their  "incompressible"  values)  was  applied  and  provided  the  same 
results.  In  Figure  15,  results  are  also  shown  with  the  J-L  model 
and  the  PTM  technique  for  cases  where  A  and  B  are  divided  by  M, 
and  where  B  alone  is  divided  by  M"^ .  In  this  case,  the  best 
agreement  is  obtained  when  A  and  B  are  divided  by  M. 

3 . 2  Sharp  cones 

3.2.1  Martellucci  experiment,  sharp  cone. 

Heat  transfer  data  on  a  7.2*  half-angle  cone  have  been 
obtained  by  Martellucci  et  al'.  (ref.  23)  . 

Figure  16  shows  comparisons  of  calculated  and  measured  Stanton 
number  for  the  following  flow  conditions: 

free  stream  Mach  number  =  7.99, 

free  stream  temperature  =  68.59*K, 

free  stream  Reynolds  number  per  foot  =  3.79x10^, 

T^T^  =  0.3.  TUE  is  set  to  1%.  The  onset  of  transition  is  imposed. 

Figure  16  displays  four  computed  curves  of  Stanton  number: 
curve  1:  J-L  model  without  PTM  technique. 

curve  2:  J-L  model  with  PTM  technique  (A  and  B  set  to  their 
incompressible  values) , 

curve  3:  J-L  model  with  PTM  technique,  A  divided  by  M,  B  divided 
by  M^, 

curve  4:  7-l  ™odel  v.'ith  PTM  technique,  A  divided  by  M,  B  divided 
by  M-^. 

Results  for  this  flow  were  inconclusive,  curve  4  shows  a 
spreading  of  the  transition  zone  but  no  peak  of  heat  transfer.  In 
the  four  cases,  agreement  between  computations  and  experimental 
measurements  is  not  very  good  and  the  above  sets  of  A  and  B  values 
for  curves  3  and  4  were  not  used  anymore. 


3.2.2  Muir  and  Trujillo  experiment,  sharp  cone. 

Muir  and  Trujillo  (ref.  24)  performed  a  series  of 
experiments  to  study  boundary  layer  transition  on  an  8  degree 
half-angle  cone.  The  experiments  took  place  in  the  U.S.  Naval 
Ordnance  Laboratory's  Hypersonic  Wind  Tunnel.  Effects  of  nose 
bluntness  were  investigated.  This  experiment  was  simulated  first 
for  the  sharp  cone  model.  The  nominal  wall  to  recovery 
temperature  ratio  is  0.6. 
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The  free  stream  turbulence  level  is  set  to  2  % .  The  onset  of 
transition  is  imposed. 

The  following  four  test  conditions  were  applied: 


cases 

M 

U  /V  Xl0~^/ft 

CD  OQ 

T^  degrees  Rankines 

1 

5 . 9 

3  . 

1092. 

2 

6. 

9.7 

1127. 

3 

6. 

17  . 

1162  . 

4 

6. 

23 . 6 

1162  . 

T^  is  the  stagnation  temperature. 

Figure  17  compares  the  results  of  Stanton  number  versus 
Reynolds  number  based  on  distance  along  the  cone  surface,  for  the 
four  free  stream  Reynolds  numbers  given  above.  The  J-L  model  is 
used  without  transition  model.  The  width  of  the  transition  zone 
is  very  narrow  in  all  cases. 

Figure  18  shows  the  same  kind  of  comparisons  when  the  J-L 
model  is  used  with  the  PTM  technique.  The  PTM  parameter  B  is 
divided  by  ,  A  is  unchanged.  The  computation  predicts  the  width 
of  the  transition  zone  very  well.  In  the  laminar  region,  the 
discrepancy  between  the  numerical  results  is  attributed  to  a 
coarse  first  cell  spacing  from  the  wall. 

Figure  19  shows  another  set  of  comparisons  between  computed 
and  measured  Stanton  number.  In  this  case,  the  J-L  model  is  used 
with  a  modified  production  term.  The  production  term  of  the  k  and 
6  equations  in  the  J-L  model  is  modified  by  the  intermittency 
factor  (formula  78,  p.  20  of  Adams  report,  ref.  25).  The 
numerical  results  predict  fairly  well  the  extent  of  the 
transition  zone,  but  the  peak  of  Stanton  number  is  smeared,  and 
for  this  reason  the  PTM  technique  was  prefered. 

Figure  20  shows  the  variation  of  Reynolds  number  based  on 
momentum  thickness  versus  the  Reynolds  number  based  on  distance 
along  the  cone  axis.  It  is  observed  that  the  onset  of  transition 
occurs  at  R  Q  =  700  for  Re/ ft  =  3,000,000  and  at  R  =  900  for 
the  other  vilues  of  Re/ft.  For  this  flow,  the  edge^Mach  number  is 
5.24.  If  we  were  to  use  a  stability  criterion  the  effect  of  TUE 
would  need  be  incorporated  in  the  criterion. 


3.3  Blunt  cones. 

Numerical  predictions  of  heat  transfer  were  performed  for 
blunt  cones  by  Martellucci  et  al.  (ref.  23)  and  Muir  and  Trujillo 
(ref.  24) .  The  numerical  simulations  of  some  of  these  experiments 
are  reported  below. 
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3.3.1  Martellucci  experiment,  blunt  cone. 

The  test  conditions  are: 

cone  half-angle  =  7.2®,  ^ 

free  stream  Reynolds  number  per  foot  =  3.79x10°, 
ratio  of  i.ose  radius  over  base  radius,  R^/Rg  =  0.02, 
free  stream  Mach  number  =  8. 

TUE  =  1%. 

Figure  21  presents  an  enlargement  of  the  nose  region  mesh. 
Figure  22  shows  typical  density  contours.  Figures  23  and  24  show 
some  entropy  contours  and  an  enlargement  of  the  entropy  contours 
in  the  nose  region. 

Figure  25  compares  the  results  of  computed  Stanton  number 
versus  streamwise  distance  with  Martellucci  et  al.  experimental 
measurements.  The  J-L  model  is  used  in  the  four  cases: 

1.  no  transition  used. 

2 .  PTM  technique  used  with  A  and  B  equal  to  their 

incompressible  values.  . 

3.  PTM  technic[ue  used  with  A  unchanged  and  B  divided  by  M^ . 

4 .  PTM  technique  used  with  A  unchanged  and  B  divided  by  M^ . 

In  Figure  25,  curve  4  seems  to  provide  the  best  predictions. 


3.3.2  Muir  and  Trujillo  experiment,  blunt  cone. 

Muir  and  Trujillo  experiment  for  a  blunt  cone  was  simulated 
for  the  four  test  conditions  mentioned  above  for  a  sharp  cone  in 
Section  3.2.2.  The  model  is  a  blunt  cone  of  nose  radius  equal  to 
0.1  inches. 

The  transition  onset  is  imposed. 

TUE  is  equal  to  2%. 

Figure  26  shows  results  of  Stanton  number  versus  Reynolds 
number  based  on  distance  along  the  cone  surface,  when  the  J-L 
model  is  used  without  the  PTM  technique.  The  width  of  the 
transition  zone  is  too  narrow. 

Figure  27  shows  results  of  Stanton  number  versus  Reynolds 
number  based  on  distance  along  the  cone  surface,  when  the  J-L 
model  is  used  with  the  PTM  technique  and  with  B  divided  by  M^,  A 
unchanged.  In  this  case,  the  agreement  between  computed  results 
and  experimental  measurements  is  very  good  with  respect  to  the 
prediction  of  the  extent  of  the  transition  zone  and  the  peak 
level  of  Stanton  number.  In  the  laminar  regime,  the  discrepancy 
between  the  computed  results  is  caused  by  a  too  coarse  first  cell 
spacing. 
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4.  CONCLUSIONS. 

A  new  transition  model  was  tested  for  hypersonic  transitional 
flows  over  flat  plates,  sharp  and  blunt  cones.  The  transition 
model  was  first  developed  for  incompressible  flows,  the  effect 
of  free  stream  turbulence  levels  is  built  in  the  model.  In  the 
applications  of  this  model  to  hypersonic  transitional  flow,  the 
model  constants  had  to  be  modified.  Attempts  were  made  to 
introduce  a  free  stream  Mach  number  effect.  In  some  instances, 
the  computed  results  agree  very  well  with  experimental 
measurements.  The  level  of  free  stream  turbulence  (or  noise)  is 
seen  to  significantly  influence  the  results.  In  the  experiments 
simulated  here,  the  free  stream  noise  level  was  not  measured 
making  difficult  to  relate  the  TUE  value  to  the  noise  level. 
However,  the  best  computed  results  were  obtained  for  the 
simulation  of  the  Muir  and  Trujillo  experiment  for  TUE=2%,  which, 
from  our  crude  conversion  analysis  of  TUE  into  noise  (or  vice 
versa) ,  corresponds  to  a  free-stream  rms  pressure  fluctuations 
divided  by  mean  static  pressure  (or  noise  level)  of  4%.  From 
Figure  4,  at  Mach  6,  the  noise  level  in  wind  tunnels  can  vary 
from  1%  to  4%,  which  is  in  satisfactory  agreement  with  the 
value  used  in  the  computations. 

The  transition  model  employs  a  stability  criterion  to 
determine  the  transition  onset.  The  Bushnell  stability  criterion 
was  found  to  perform  approximately  well,  but  this  criterion  needs 
to  be  refined.  It  is  necessary  to  introduce  in  the  stability 
criterion  the  effect  of  free  stream  turbulence  level  (or  noise) . 

An  additional  contract  would  be  required  to  pursue  further 
investigations  on  this  important  aspect  of  transition  modeling. 

The  reported  results  are  encouraging  and  the  PTM  transition 
model  shows  potential  of  becoming  a  predictive  tool  for 
hypersonic  transitional  flows.  More  work  is  required  to  compare 
computations  and  experimental  measurements  for  a  large  variety  of 
hypersonic  transitional  flows  over  cones,  which  will  establish 
the  universality  of  the  model. 
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Velocity  fluctuaiions,  u 

Acoustic  Noise.  p 

Temperature  fluctuations,  T 
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Transonic 
a  6  <  M(J,  <  L  3 

Velocity  fluctuations.  u 
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Temperature  Fluctuations.  I 
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Figure  1.  Flow  disturbances  in  Wind  Tunnels. 
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Figure  2.  Variation  of  "A"  v'-'t..  free-stream  turbulence 
intensity  for  the  Jones-Launder  model. 
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Figure  3.  Variation  of  'B"  with  free-stream  turbulence 
intensity  for  the  Jones-Launder  model. 
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Free— stream  rms  pressure  fluctuations  divided  by 
mean  static  pressure  versus  Mach  number. 
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Figure  5.  Skin  friction  versus  Reynolds  number  based  on  streamwise 
distance,  computed  by  Schmidt  and  Patankar  with  the 
Jones-Launder  model  without  the  PTM  technique. 


Figure  6.  TURF  computed  skin  friction  versus  Reynolds  number 

based  on  streamwise  distance,  with  the  Jones-Launder 
model  and  without  the  PTM  technique. 
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Figure  7.  Skin  friction  results  of  Schmidt  and  Patankar 

obtained  with  the  Jones-Launder  model  and  the  PTM 


technique . 
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Figure  8.  TURF  computed  skin  friction  with  the  Jones-Launder 
model  and  the  PTM  technique. 
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Figure  13.  Comparisons  between  computed  and  experimental 
heat  transfer,  R  /in.  =  560,000. 
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Figure  15.  Comparisons  between  computed  and  experimental 
heat  transfer. 
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Figure  17.  Streamwise  variation  of  Stanton  number,  no 
transition  model  used. 
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Figure  18.  Comparisons  between  computed  Stanton  number  and 
experimental  measurements.  The  PTM  technique  is 
used. 
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igure  19.  Streamwise  variation  of  Stanton  number,  the 

intermittency  factor  approach  of  Adams  is  used 
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Figure  20. 


Variation  of  Reynolds  number  based  on  momentum 
thickness  versus  R 
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Figure  23.  Entropy  contours,  blun':  cone  at  Mach  8. 
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Figure  24.  Enlargement  of  entropy  contours,  blunt  cone  at 
Mach  8 . 
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Stanton  number 


Figure  25.  Calibration  of  A  and  B  in  the  PTM  model. 
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Figure  27.  Variation  of  Stanton  number  with  free-stream  Reynolds 
number  for  a  blunt  cone,  RN=0.1  in.  -  PTM  model  used. 
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